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Abstract 

We calculate the one-loop corrections to a general off-forward deeply-virtual 
Compton process at leading twist for both parton helicity-dependent and in- 
dependent cases. We show that the infrared divergences can be factorized en- 
tirely into off-forward parton distributions, even when one of the two photons 
is onshell. We argue that this property persists to all orders in perturbation 
theory. We obtain the next-to-leading order Wilson coefficients for the general 
leading-twist expansion of the product of two electromagnetic currents in the 
MS scheme. 
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I. INTRODUCTION 



Photons, real or virtual, are known to be clean probes of the internal structure of the 
nucleon. In deep inelastic scattering (DIS), the cross sections for absorption of highly- virtual 
photons were the first to reveal the internal quark structure of nucleons. The parton dis- 
tributions extracted from these cross sections contain important structural information and 
seriously challenge our understanding nonperturbative quantum chromodynamics (QCD). 
Elastic absorption of virtual photons can be used to measure the electromagnetic form fac- 
tors of the nucleon. At low virtuality, these form factors give us direct information about 
the sizes and magnetic moments of nucleons. At high virtuality, they are sensitive to the 
leading-twist light-cone wavefunctions. More recently, real photon elastic scattering at low 
energy has been used to extract the electromagnetic polarizablities of nucleons. 

In a recent paper, one of us introduced deeply-virtual Compton scattering (DVCS) as a 
probe to a novel class of "off-forward" parton distributions (OFPD's) ][]]. DVCS is a process 
in which a highly virtual photon (with virtuality Q 2 >>> Aq CD ) scatters on a nucleon target 
(polarized or unpolarized), producing an exclusive final state consisting of a high-energy 
real photon and a slightly recoiled nucleon. With the virtual photon in the Bjorken limit, a 
QCD analysis shows that the scattering is dominated by the simple mechanism in which a 
quark (antiquark) in the initial nucleon absorbs the virtual photon, immediately radiates a 
real one, and falls back to form the recoiled nucleon. 

Several interesting theoretical papers have since appeared in the literature, which studied 
the DVCS process further. In Ref. @, the single-quark scattering was recalculated using a 
different, but equivalent definition of the parton distributions. The evolution equations of 
the distributions were derived and some general aspects of factorization were discussed. In 
Ref. ||, the evolution equations for OFPD's were derived and the leading- twist DVCS cross 
sections were calculated at order Some past and recent studies of OFPD's can be found 
in ||. In Ref. £§], estimates of these cross sections were made at COMPASS and TJNAF 
energies. In Ref. [f| , the DVCS process was considered as a limit of unequal mass Compton 
scattering, which was studied from the point of view of the operator product expansion. 
Some early studies of unequal mass Compton processes can be found in Refs. |7]||. In Ref. 
||, a number of suggestions were made to test the leading twist dominance in DVCS at 
finite Q 2 . In a Rapid Communication paper, the present authors studied 0(a s ) corrections 
to DVCS for the parton helicity-independent case QT0j| . In Refs. 0,0, the same issue was 
investigated from different perspectives. The present paper is an expanded presentation of 
our results in Ref. JTIJ . 

The main motivation for the present study is to see if the theoretical basis for the DVCS 
process is up to par with other well-known perturbative QCD processes. More explicitly, 
we discuss the existence of a factorization theorem for this process. For general two virtual 
photon processes in the Bjorken limit, the factorizability is suggested by studies of deep 
inelastic scattering. In the case of DVCS, where one of the photons is onshell, the situation 
could be different. Potential infrared problems can arise because of the additional light-like 
vector in this special kinematic limit. However, it is believed that these complications will 
not ruin the factorization properties |§. 

To see factorization at work, it is instructive to work out one-loop examples. We will do 
this explicitly in section III. For consistency, we consider the unphysical process of DVCS on 
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onshell quark and gluon "targets." To ensure gauge invariance, we regularize the infrared 
divergences by going to d = 4 + e dimensions. For completeness we have considered both 
the symmetric and antisymmetric parts of the amplitudes, which are related to helicity- 
independent and dependent parton distributions, respectively. The only omission is the 
gluon helicity flip amplitude, which will be discussed in Ref. As expected, our result 
contains collinear infrared divergences which can be interpreted as the one-loop perturbative 
parton distributions, as we will show in Section IV. This property is independent of the 
special kinematic limit of DVCS. 

A general proof of the DVCS factorization was first given by Radyushkin in his approach 
based on a-representation ||. In this paper, we give an alternative proof using the tools 
developed by Libby, Sterman, Collins and others Hl4f . According to these, one can represent 
the infrared sensitive contributions in a generic Feynman diagram with reduced diagrams. 
These reduced diagrams have intuitive physical significance and are easy to identify. General 
power counting rules can be used to select leading reduced diagrams in a process. A recent 



application of the method can be found in Ref. [15] . We show in Section V that the leading 
reduced diagrams for DVCS do not contain any soft divergences and are in fact exactly 
the same as those present when the final state photon is deeply virtual. The collinear 
divergences in the reduced diagrams can be attributed to those of OFPD's when calculated 
in perturbation theory. Therefore we conclude that factorization for DVCS is in the same 
footing as other well-known examples like deep inelastic scattering. 

The factorization property of the general two virtual photon process can be summarized 
beautifully in terms of Wilson's operator product expansion. This expansion requires op- 
erators with total derivatives [H-H] to describe the off-forward nature of the process. It is 



well-known that these derivative operators contribute to the wavefunctions of mesons ||16|| . In 
section VI, we convert our one-loop results into Wilson coefficients of the twist-two operators 
in the MS scheme. Together with the two-loop anomolous dimensions of these operators, 
they provide the necessary ingredients for calculating DVCS at the next-to-leading order. 
We summarize and discuss our results in Section VII. 



II. KINEMATICS AND PARTON DISTRIBUTIONS 

Although our ultimate interest is in deeply virtual Compton scattering, we start by con- 
sidering a general Compton process involving two offshell photons with different virtualities. 
This and a suitable choice of kinematic variables allows us to exploit the full symmetry of 
the problem. In the general Compton process, a virtual photon of momentum q + A/2 is 
absorbed by a hadron of momentum P — A/2, which then emitts a virtual photon with 
momentum q — A/2 and recoils with momentum P + A/2. The three independent external 
momenta can be expanded in terms of the light-cone vectors 

p" = (p + ,0,0,p + ) , (1) 
^ = ^(1,0,0,-1) , 

where the 3-direction is chosen as the direction of the average hadron momentum (P), and 
two transverse vectors. In an expansion, we call the coefficient of p M the + component and 
that of the — component. Thus we write 
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M 2 — tlA 
2 

<f = -CP" + , (2) 

A M = -2&" + £(M 2 - f/4)ra" + A£ , 

where M is the hadron mass (which is taken to be the same for the initial and final hadrons), 
t = A 2 , Q 2 is the virtuality of g M , £ is a measure of the difference of the virtualities of the 
two external photons, ( is defined as 



Q 2 



-1 + 



1+ 4*|(M*- t /4) 



s 2i B (M2 - i/4) ^ M Q 

and Aj. is a vector in the transverse directions which has squared length —t (1 — £ 2 ) — 4£ 2 M 2 . 
We have also introduced xb = Q 2 /(2P ■ q), the analogue of the Bjorken scaling variable in 
this off-forward process. We note that these expressions limit the range of £ to 

* ~ -t + 4M 2 ' y } 



for fixed t, or the range of t to 



4£ 2 M 2 

~~ 1 - l^t2~ ' ( 5 ) 



for fixed £. 

In the Bjorken limit, these expressions simplify considerably. Since we consider only the 
leading twist in this paper, we may neglect all but the + components of P M and A M (in order 
to form large scalars, one must dot the + component of a vector with the — component of 
q). Hence, in the limit Q 2 — > oo (t remaining finite), we may write 

<f ~ -XbI? + , (6) 

ZX R 



-2ft/ 1 . 



Here, we note that the external invariants have been reduced from six to three by enforcing 
kinematics and taking the Bjorken limit. We express these three scalars in terms of one 
mass scale, Q 2 , and two dimensionless parameters, xb and £. When we introduce the parton 
distributions, our expressions will also involve the parton light-cone momentum variable x. 
Hence, the final result will be expressed as algebraic functions of xb, x, and £ multiplied by 
the appropriate power of Q 2 . 
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FIGURES 
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FIG. 1. The handbag diagram for the general two photon process. 



Our goal is to factorize the short and long distance physics of the Compton amplitude 
in the Bjorken limit, 
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d ze 
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(7) 



where J M = J2 q Zqtyql^tyq is the electromagnetic current and ip q is the bare quark field of 
flavor q and charge e q . The simplest Feynman diagram for this process is shown in Fig. 1, 
where a quark comes out of the nucleon blob, scatters, and rejoins the nucleon blob. While 
the scattering involves a large momentum transfer and can be calculated in perturbation 
theory, the nucleon blob with two quark legs is related to the baryon structure and is 
nonperturbative. For more complicated graphs, as will be discussed throughout this paper, 
the Compton amplitude can be separated analogously into soft and hard contributions. In 
the remainder of this section, we will highlight some important aspects of the soft part. 

The nonperturbative contribution to the Compton amplitude in Eq. (0) can be expressed 
in terms of off- forward parton distributions contained in the parton density matrices ||. For 
quarks, we call the density matrix M a p, where a and (3 are Dirac indices, and expand it in 
terms of the Dirac matrices. At leading twist, M a p is just the light-cone correlation function, 



2vr \ 2 



\ F q ^ a p + \F q ( 75 + 



-iq f Y , n-A(Cn)dC 





(8) 



where the ellipses denote contributions either of higher twist or chiral-odd structure, which 
do not contribute to the leading process under consideration. The P symbol denotes the path 
ordering of the exponential, which makes this expression gauge invariant. It is necessary 
to include this gauge link whenever one is not working in the light-cone gauge (A + = 0). 
Multiplying by fip a and (7^75)^ and taking traces, we project out the same distributions 
as considered in [Hi: 



Fq(x,0 
Fq(x,0 



1 rdX 
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(9) 
(10) 
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We have suppressed the renormalization scale \i which is always present in defining a parton 
distribution. We have also suppressed the t dependence because it will not affect most of 
the discussions in this paper. 

At next to leading order, gluons also contribute to the Compton process. Although 
it is nontrivial to show, the twist-two gluon distributions are contained in the following 



gauge-invariant light-cone correlations (e 
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(11) 



c »a~ ( 6 p v p f3 _ e ^yS p u p a + e ^jS pfip a _ ^ayS^^ ^ + 



= -xF G (x,0 (g^f^ 
+ixF G (x,£) (e 

where the ellipses denote higher twist contributions and an additional twist-two term which 
involves gluon helicity flip and will not be considered in this paper [ 13| . Again, P{} denotes 
path ordering (we note that here the gauge link is in the adjoint representation of SU(3)). 
The off-forward gluon distribution functions Fq and Fq may be isolated by contraction and 
are 
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(12) 
(13) 



Here, we have defined the dual field strength tensor F^ u = \e^ va ^F t 



a(3- 



It is easiest to see the connection of the above gluon distributions with the nonpertur- 
bative structure arising from Feynman diagrams in the light-cone gauge. In this gauge, the 
gauge link is just the unit operator in the adjoint representation and field strength tensors 
with one + index F +fl simplify to d + A fJ '. Fourier transformation to momentum space yields 
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(14) 



(15) 



where we have defined x + = x + £ and x_ = x — £. VT represents (2ir) 5 (0), the space- 
time volume of our system. In a factorized calculation of the Compton amplitude involving 
gluons, the gluonic indices in the hard part will be contracted with the tensor 

d 4 < 



_L / JlL x / P + ^ 

VTi f27r) 4 V ; \ 2 
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(16) 
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In the above definitions, we have assumed that we are working in 3+1 space-time di- 
mensions. However, to regularize the ultraviolet and infrared divergences arising from loop 
diagrams, it is convenient to generalize them to d dimensions. Let us first consider the 
quark density matrix in Eq. (8). Because the spinors are kept in 4 dimensions, the first 
term on the right hand side generalizes to d dimensions without change. The second term, 
however, involves 75 which has no unique extension. Different choices, in the end, define 
different factorization schemes. If one uses the t' Hooft-Veltman definition (75 = i / -f 0/ -f l/ -f 2/ y 3 ) 



17 1, one usually introduces an extra renormalization constant Z 5 so that the non-singlet 
axial currents are conserved. An alternative choice is offered by Chanowitz, Furman, and 
Hinchcliffe ||18|| , which employs the usual four- dimensional rules 

{7^,751 = tyeM 
Tr[ 757 V7V] = -(4 + 0(e)) ie a ^ 8 . (17) 

The ambiguity in the second equation does not affect calculations as long as there are no 
anomalies in the problem. In the case that there is an anomaly, the ambiguity can be fixed 
by imposing the relevant Ward identities. 

We now turn to the gluon density matrix in Eq. (TO). Fq{x) contains an average over 
gluon polarizations. To make this consistent with the number of transverse polarization 
states available to gluons in d dimensions, we multiply this term by 1/(1 + e/2). The 
polarized gluon density is related to the antisymmetric combination of the gluon fields F +l 
and F +2 . This does not change after going to d dimensions if the target polarization is kept 
the same. Hence, we have left that term as it is. 



III. ONE-LOOP COMPTON AMPLITUDES 
ON QUARK AND GLUON "TARGETS" 

In this section, we present a one-loop calculation of the general Compton scattering on 
onshell quark and gluon "targets" in the Bjorken limit. The result will be used in the next 
section to show that the factorization of soft and hard contributions can be done consistently 
with the definition of off-forward parton distributions. In particular, this property does not 
change in the limit of a real final state photon. Our result will also be used to derive a 
generalized operator product expansion to next-to-leading order. For the convenience of 
the reader, we are going to spell out some technical details of the one loop calculation. 
We believe that some of the techniques, like the cancellation of propagators and light-front 
coordinate integration, will be useful in other contexts. 

We begin with an onshell quark "target". Here, there are two diagrams at leading 
order (LO) and eight at next-to- leading order (NLO). Half of these diagrams are shown in 
Fig. 2. The other half will be taken into account by using the crossing symmetry, i.e., 
the simultaneous replacement of q — > — q and //•«-> v. The terms with Xb — > — %b i n 
the following formulas reflect this contribution. Because of time reversal invariance, the 
Compton amplitude is also an even function of £, i.e., symmetric under £ — > — £. This 
symmetry relates the left and right vertex diagrams (Figs. 2c and 2d, respectively) to each 
other. On the other hand, the quark self-energy diagram and the box diagram are themselves 
^-symmetric. This symmetry not only allows us to reduce the number of graphs at NLO 
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from four to three, but also becomes a powerful tool which helps us compute each amplitude, 
as we illustrate later. 




(b) (c) (d) (e) 

FIG. 2. Diagrams for Compton scattering on an onshell quark to order a s . 



Before presenting the details of our calculation, it is necessary to discuss the issues of 
onshell reduction and ultraviolet divergences. We calculate diagrams in d dimensions and use 
Feynman's gauge. As such, we can take the quark target to be massless and take the onshell 
limit at the beginning of our calculation. In the modified minimal subtraction scheme (MS), 
the renormalized single-quark propagator has a residue, SZ, at the pole p 1 = 0. According to 
Lehmann-Symanzik-Zimmerman reduction formula, we calculate the onshell physical matrix 
element of an operator O r as 

(p\6 R \p) = \im 5ZJ2G , (18) 

where Go is the set of all amputated connected graphs with one insertion of Or in renor- 
malized perturbation theory. The factor 5Z is infrared divergent for massless quarks and 
equals Zp 1 in the present calculational scheme. Since J M = Zp^jjR'j^ipR, all renormaliza- 
tion constants, including the subtraction for the quark self-energy, cancel at one-loop level. 
Therefore, T^ v for single quark and gluon "targets" can be calculated just from the graphs 
shown in Fig. 2. 

Examining these graphs, we see that the self-energy diagram contains a loop integral 
with two Feynman denominators, the vertex diagram contains one with three, and the 
box with four. A one-loop integral with two propagators is straightforward. Difficulties 
arise, however, with the calculation of three and especially four-propagator integrals. These 
difficulties may be avoided in this calculation because of several simplifications. Consider 
first the box diagram. The loop integral is of the form (with momentum routing as shown 
in Fig. 3) 

r d d k Tr[^^-^)Y^+ij)r^ + U)la^ 5 )} nq) 
J (2ir) d (k + ip) 2 {k - ip) 2 {k - pf{k + q) 2 ' 1 ' 

where we have replaced A M by — 2£p M . In order to simplify the integral, we express the trace 
as a sum of terms which cancel one of the propagators. This can be done because both k 2 
and 2p ■ k can be written as linear combinations of (k + ^p) 2 and (k — £p) 2 , and the trace 
vanishes whenever k~ and k\ do. 
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k-?p 



(1-S)P 



FIG. 3. The momentum flow in the box diagram 



Now that we have shown that a denominator can be cancelled, we have effectively re- 
duced the four-propagator problem to a three-propagator one. Since we have only shown 
that the numerator will be a linear combination of two different denominators, rather that 
proportional to one, the four-propagator integral will in general become two three-propagator 
integrals. However, we may use the £ symmetry by writing 



k 2 

2p ■ k = 



(20) 



In this way, we consider only the (k + £p) 2 cancellation and let the symmetry take care 
of the rest. We note that the integral we now have is exactly the same (up to numerator 
differences) as that arising from the right vertex correction. We will see that this basic 
integral is the only one we must calculate to obtain both the polarized and unpolarized 
amplitudes for both the quark and gluon contributions (if we forget, for the moment, the 
simple self-energy diagram). The integral has the form 



d d k 



Numerator 



(2n) d [(k - p) 2 + ie] [(k - tpf + ie] [(k + q) 2 + ie] 



(21) 



We have found that this integral is easily done in light-cone coordinates by expanding k^ in 
terms of the light cone momenta, p^ and n M . We first do the k~ integration by contour in 
the unphysical region of large xb, and then the transverse integrations. The k + integration 
is left until the end. If we write k + = yp + , the value of the integral fl2l|) is 
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where we have defined N according to 
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a + (3k 



N = (3 — 
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a xb(xb — cl) 
Q 2 (y ~ a){x B - y) 



(22) 



(23) 



Doing the y-integrals requires some care because a delicate cancellation must occur if one is 
to get finite result, but the treatment is straightforward. After the £ and crossing symmetries 
are used, we find the full NLO result for the symmetric quark amplitude 
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In 1 - 



Or 



(24) 



x B — 1 

and the antisymmetric amplitude 



9' 




- - (25) 



x B 



X B 



Here we have introduced Cp = Tr[t a t a ] = N ^~ l in SU(N C ), where N c is the number of 
colors. We note that the divergences in these amplitudes are, in fact, infrared divergences 
since renormalization has removed all ultraviolet singularities. Their presence signals the 
existence of nonperturbative physics in the process. As mentioned earlier, these divergences 
will be factorized into nonperturbative matrix elements whose values can be extracted from 
experiment. We will explicitly show this in the next section. For now, we summarize the 
results of the gluon piece of the calculation. 



q-5p 



d-5)P 





(l+l)P 



(a) (b) 

FIG. 4. Diagrams for gluon Compton scattering at one loop. 

There are six graphs which contribute to the LO amplitude for gluon-photon scattering. 
These six can be reduced to three by reversing the fermion number flow, and one of these 
can be eliminated by crossing symmetry. The two distinct graphs we must calculate are 
shown in Fig. 4. The denominator of Fig. 4a is identical to that of the quark box. Again, 
the numerator is seen to vanish whenever k 2 and 2p ■ k do, which allows us to cancel one of 
the propagators exactly as above. Fig. 4b is somewhat more tricky. This diagram is itself 
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symmetric under both crossing and £ symmetry. Labeling the momenta as shown, we see 
that under an integral the symmetry q — > — q is equivalent to A <-> B and C <-> -D and 
£ — > — £ is equivalent to B ^ D and A <-> C '. We also note that here it is possible to 
represent 1 as a linear combination of Feynman denominators, which guarantees our ability 
to cancel one. Of course, since 1 does not depend on q or £, it may be represented in the 
symmetric form 



1 = 



2(1 -Op -g 



A 2 + (£--0 + (?->-?)■ 



(26) 



Now it remains only to calculate the trace in a symmetric way and substitute the result into 
the formulae of the quark calculation. Averaging the gluon polarization for the symmetric 
amplitude in (i-dimensions, one finds 
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for the symmetric amplitude and 



^ 27T 



g 2 e 



i 



- 2x B 

- {x B 



i-e \ 



i + 4 

e 



In 



x B 



(28) 



for the antisymmetric amplitude. We have defined Tp = Tr[t a tb] = |5 a f,. 



IV. ONE-LOOP FACTORIZATION 
AND EVOLUTION OF OFF-FORWARD PARTON DISTRIBUTIONS 

We now turn to the infrared divergences present in all of the amplitudes. These di- 
vergences arise from the regions of loop-mementum integration where some of the internal 
propagators are near their "mass shells". In these regions, perturbative calculations are 
clearly meaningless. The standard procedure of fixing this problem is to factorize the ampli- 
tudes into the infrared safe (i.e. devoid of infrared divergences) and infrared divergent pieces, 
interpreting the latter as nonperturbative QCD quanities. Of course, the factorization pro- 
cedure has a large degree of arbitrariness. To factorize in a physically interesting way, one 
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usually chooses the nonperturbative objects as the parton distributions in the target, defined 
in a particular renormalization scheme. In this paper, we consider parton distributions in 
the MS scheme. The goal of this section is to show that all infrared divergences present in 
the Compton amplitudes can be associated with these distributions. 

Since we consider factorization within the framework of perturbation theory, we need to 
compute the parton distributions in quark and gluon targets in perturbative QCD. At the 
leading order in a s , one has 



F° /q (x,0 = 6 qq ,6(x-l) , 



(29) 



in an onshell quark q. At next-to-leading order, one can calculate directly from the definitions 
in section II (with the external hadron states replaced with perturbative quark states): 
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The quark distributions contain infrared divergences signaled by the presence of the 1/e 
terms. These divergences reflect the soft physics intrinsic to the parton distributions. F^ g , 
calculated for the quark target, satisfies the evolution equation derived in Ref. : 
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where b is summed over all parton species and the P's are the off-forward Alterelli-Parisi 
kernels, or splitting functions. The 'covariant' derivative is defined to include A, the endpoint 
contribution in Eq. fl3~0|). 

We can reexpress the symmetric part of the quark Compton amplitude in terms of the 
unpolarized, off-forward quark distribution 
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Analogously, we find that the antisymmetric part of the quark Compton amplitude can be 
expressed in terms of the polarized off-forward quark distribution 
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We now turn to Compton scattering on a gluon "target" . Infrared divergent contributions 
come from the regions where the quarks in the box diagrams are nearly onshell. Therefore, 
it is natural to associate these divergences with the quark distributions in a gluon target. 
On the other hand, the finite contributions come from regions where large momenta run 
through the quark loop. In these regions, the photon has an effective pointlike coupling 
with the gluons in the target. At leading order, the off- forward gluon distribution in a gluon 
target is just 



x ,Z) = -(5(x-l)-5(x + l)) . 



(36) 



To order a s , there are quark partons in the gluon. The corresponding off-forward quark 
distribution is 



a s (fx) (2 



2tt 



- + ln(^e'»/47r^))P(x,0 , 



where for x > £ 



P(x,£) = 27> 



x 2 + (l 



X) 



e 



and for — £ < x < £ 



P(x,0 = T F - 



(i-e) 2 

x + £)(l-2:r + £) 



£(i+o(i -a 



(37) 



(38) 



(39) 



F} /g (x,0 = Oforx<-£. 

With the above off-forward distributions, we reexpress the symmetric part of the gluon 
Compton scattering amplitude 
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(40) 
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Similarly, we can reexpress the antisymmetric part of the gluon Compton amplitude in terms 
of helicity-dependent, off-forward quark and gluon distributions 



9 



~ F a3ii 

i— — e pj n a pf3 




(41) 



Xb — x 



We summarize the one-loop Compton amplitude on a target N in the factorization 
formula 



2 N 
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1 dx 
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-)\ 


\x B 





(42) 



where C q , C g , C q and C g are shown to order a s in Eqs. (H H, ||, 0), respectively. [We 
emphasize here again that we have neglected the contributions from longitudinally polarized 
photons and from photon helicity flip. Both effect start at order a s .} In the above form, 
all infrared sensitive contributions have been isolated in the relevant parton distributions, 
which must be calculated nonperturbatively or measured in experiments. In the DVCS limit 
£ — > xb, the coefficient functions remain finite, although they have branch cuts there. This 
indicates that factorization holds for two-photon amplitudes even when one of the photons 
is onshell. We will argue in the next section that the above formula, one of the main results 
of this paper, remains valid to all orders in perturbation theory. 



V. FACTORIZATION OF DVCS AMPLITUDES TO ALL ORDERS 

In this section, we generalize the one-loop result of the previous section, showing that 
the factorization formula Eq. fl4*2| ) is valid in the DVCS limit to all orders in perturbation 
theory. The one-loop result indicates that all soft divergences - those associated with inte- 
gration regions where all components of some internal momenta are zero - cancel, whereas 
all collinear divergences can be factorized into the off-forward parton distributions. To see 
that this happens also at higher orders in perturbation theory, it is important to understand 
how the soft cancellation happens in the simplest case. 

The self-energy diagram in Fig. 2b does not contain any infrared divergences because 
the intermediate quark is far offshell. The vertex corrections in Figs. 2c and 2d potentially 
have infrared divergences, but a simple power counting indicates that these diagrams are in 
fact infrared convergent. Thus infared divergences appear only in the box diagram. In the 
region where the gluon momentum k is soft, we can approximate the integral as 
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d*k , 111 

llf'Pzr-L—Lui* ( 43 ) 



(2vr) 4 ^ p'-kp-kk 2 

where p and p' are the momenta of the two external quark lines. On the other hand, the 
wave function renormalization of an "on-shell" quark, SZ, also contains infrared divergences. 
Grouping these divergent terms together, we have the entire soft contribution 



1 r d A k ( p» p'» V 1 



2 J (2vr) 4 \p ■ k p 1 -k 



(44) 



In the collinear approximation, p' is proportional to p and thus the above integral vanishes. 

For higher-order Feynman diagrams, a systematic method of identifying, regrouping and 
factorizing infrared-sensitive contributions has been developed by Libby, Sterman, Collins, 
and others [TJ]]. The method essentially consists of the following steps: 1) simplify the 



Feynman integrals by setting all the soft scales to zero, including the quark masses; 2) 
identify the regions of loop integration which give rise to infrared divergences; 3) use infrared 
power counting to find the leading infrared-divergent regions; 4) show that all soft and 
collinear divergences either cancel or factorize into some nonperturbative quantities. In the 
remainder of this section, we examine the validity of the factorization formula Eq. (^) in 
the limit of xb = £ following the above steps. 

In DVCS, the leading contributions come essentially from a massless collinear process in 
which the external momenta take the form shown in Eq. (6). In this simplified kinematic 
region, all infrared-sensitive contributions appear as 1/e poles in dimensional regularization. 
If a contribution contains no infrared divergences, it comes from regions of loop momenta 
comparable to the hard scale Q 2 , and thus it is insensitive to the soft scales. An infrared 
divergent contribution must come from the integration regions where some internal propa- 
gators are near their mass shells. Since such soft contributions cannot be calculated reliably 
in perturbative QCD and eventually must be taken into account with nonperturbative ma- 
trix elements, one can use any valid infrared regulator to characterize them in perturbative 
calculations. Thus the collinear massless limit helps to simplify the identification of soft 
contributions while leaving the truly-perturbative contributions intact. 

Infrared divergences appear in a Feynman diagram when some of the external momenta 
are onshell. The regions of integration producing such contributions can be identified from 
the Landau equations which are derived by considering the analytical properties of the 
diagrams as functions of complex external momenta. According to Coleman and Norton 
l9| , these regions can be represented by the so-called reduced diagrams in which offshell 



lines are shrunk to points and onshell lines are drawn according to their real space-time 
propagation. We shall argue below that the general leading reduced diagram for DVCS is 
the one shown in Fig. 5a, in which an incoming virtual photon and an outgoing real one 
are attached to the hard interaction blob, which in turn is connected to the forward nucleon 
jet with two collinear quark lines or two physically polarized gluon lines, plus an arbitrary 
number of longitudinally-polarized collinear gluon lines. 
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(a) (b) (c) 

FIG. 5. General reduced diagrams for the DVCS process. 



To decide that a reduced diagram is leading, one can use infrared power counting J20 
which is essentially a light-cone dimensional analysis [IB]. A simple way to proceed is to 
consider the mass dimensions of the soft vertices that connect lines with either collinear or 
soft momenta. Since the dimension of an amplitude is fixed, all soft mass dimensions must 
be compensated by the hard scale Q 2 . Assuming covariant normalization for the external 
states {(p\p) = 2p°(2ir) 3 S 3 (0)), every external wave function contributes mass dimension — 1. 
The collinear quarks and gluons into a soft hadron vertex have effective mass dimensions 
depending on their polarizations. A Dirac field ip can be written as a sum of good (t/>+) 
and bad (■?/>_) components, where ip± = P±ip and P± = |7 =F 7 ± . The good (bad) component 
has effective light-cone mass dimension 1 (2). A vector potential A M has light-cone com- 
ponents A + , A , and A~ , which have effective mass dimensions 0, 1, and 2, respectively. 
For the reduced diagram shown in Fig. 5a, the only soft mass dimension comes from the 
nucleon-quark-gluon blob. Using the above rule, we find it is [0 = 2 (physical parton lines 
— 2( external nucleon states)]. Because T^ u is dimensionless, the leading reduced diagram 
contributes at order O(Q ). 

It is somewhat surprising that the leading region is independent of the virtuality of the 
final state photon as long as the initial photon is deeply virtual. When the final state photon 
is real, it can have pointlike coupling to quarks as well as extended coupling via its soft wave 
function, as happens in the case of vector dominance. Thus one has an additional reduced 
diagram in which a jet of quarks and gluons emerges in the direction and combines into 
a real photon long after the hard scattering (Fig 5b). Such a reduced diagram has already 
been considered in Ref. [1~5"| and is 0{Q^ 1 ) by infrared power counting. Indeed, according 
to the discussion in the previous paragraph, the photon wave function vertex has a soft 
mass dimension 1 = 2 (quark lines) —1 (photon state) (recall the dimensionful pion decay 
constant f n ). A negative hard power (Q _1 ) is needed in T^ u to balance it out. 

The power counting involving soft quark and gluon lines is more subtle, and some dis- 
cussion may be found in Ref. |15|]. The result is that any reduced diagram with soft lines 
connecting the hard scattering blob to the nucleon jet is subleading (Fig. 5c). The situa- 
tion here is exactly analogous to the case of forward virtual Compton scattering relevant to 
deep- inelastic scattering as discussed, for instance, by Sterman p0fl . A simple example is 
the vertex correction diagram we already discussed above. When the gluon becomes soft, it 
has a reduced diagram like Fig. 5c. However, it has no contribution at leading power. 

We now come back to the leading reduced diagram shown in Fig. 5a. The collinear 
gluons with longitudinal polarization connect the hard scattering part and the nucleon jets. 
These collinear gluons can be factorized using the generalized Ward identities as described 
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in Ref. ||14|| . Eventually, all collinear gluons can effectively be attached to an eikonal line 
in the light-cone direction conjuate to the nucleon jet, n M . Physically, this means that as 
far as the collinear gluons are concerned, the hard interaction part acts as a jet of particles 
propagating along n M . The internal structure of the hard interaction cannot be resolved 
and thus only the total color charge and momentum of the jet is relevant. The eikonal 
line together with the physical quarks and gluons and the nucleon jets form the off-forward 
parton distributions defined in Section II. In the hard scattering, only the total momentum 
and charge supplied by collinear partons are important. Thus, one can calculate it with 
incoming physical partons carrying the total momentum and charge of all the collinear 
longitudinally-polarized gluons. In this way, we have a complete factorization of the soft 
and hard physics in the DVCS process. 



VI. GENERALIZED OPERATOR PRODUCT EXPANSION 
AND WILSON'S COEFFICIENTS TO THE NLO ORDER 



The factorization formula in Eq. fl4"2| ) for the general Compton scattering process can also 
be examined in the form of an operator product expansion. The OPE was first introduced 
by Wilson [^TJ in 1969 and has been used extensively in deep inelastic scattering and other 
perturbative QCD processes. For the product of two currents separated near the light-cone, 
the expansion is threefold. Primarily, it is a twist expansion, in which twist-two contributions 
are leading whereas the higher twist terms are suppressed by powers of 1/Q 2 - Each term in 
the twist expansion contains an infinite number of local operators of the relevant twist. This 
may be thought of as a kind of Taylor expansion of bilocal operators along the light-cone. 
Finally, the coefficients of local operators (Wilson coefficients) are themselves expansions in 
the strong coupling constant. The Wilson coefficients for the unpolarized DIS process were 
calculated at order a s in the MS scheme in |22]. For the polarized case, one can find them 



m 



23|. 



When considering off-forward processes, the expansion of an operator product must 
include operators with total derivatives. We call this expansion the generalized OPE. In the 
remainder of this section, we will recast our factorization formula in its generalized OPE 
form. In the process, we identify these total derivative operators and obtain their Wilson 
coefficients to next-to-leading order in a s . The final result agrees with the known OPE in 
the DIS limit (f 0). 

To derive the generalized OPE, we expand T 4 - 7 in a power series about xb = oo. In 
this way, we can express the amplitude in terms of moments of the parton disributions 
rather than the distributions themselves. Eventually, we will relate the moments of parton 
distributions to the matrix elements of local operators. After the aforementioned expansion, 
we have 



T' 



neven=2 rneven^^" 1 X \XbJ \%B 

t a = z E E / — — 

nodd=l meven=0 J_1 X ^ Xb 




^limFqix, £) + (? nm Fg(x, 



(45) 
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The coefficients for the moments of the quark distributions in the expansion are 



S m n — 



An 



9-- + 

n n + 1 



-S 1 (n-1)[3 + - + -^— 
n n + 1 



+ AS 2 (n-l)-AT?(n-l) 

Qn 1 



>m0 



—SAm - 1) - 251(71 + m-l) 
to 



m(m + n) (n + m)(n + to + 1) 

1 1 - -)}, 



n + m n + m + 1 



'1-5, 



mO) 



(46) 



<5m0 — 



a s C F 
An 



6 



9 - - + AS 2 (n - 1) - ATl(n - 1) 
n 



-S 1 (n- 1) (3 + - + 



-251 (n + m- 1) 
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(47) 



(48) 



Notice that the above expansion contains only positive powers of x and £. This result is not 
immediately obvious because of the x + X- denominators in the amplitudes. In the case of 
the gluon distribution functions, we have an additional factor of x + X- in the denominator. 
Since the final OPE contains only local operators, these factors have to be cancelled in the 
process of expansion. Indeed, this turns out to be the case and we obtain the coefficients of 
the positive moments of the gluon distributions 
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^ m = ^ £ [2 + 25 1 (n + m-l)]( 
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Having obtained an expansion involving the moments of the distributions, we move 
toward a general form of the OPE. To this end, we consider the moments of the parton 
distributions. We begin by observing that for quarks 



1 i dxx n - 1 F q (x^) = \{Pf 



Vv(0) id^ ■■■id fMn _ 1 7^(0) 



Pi ) n 



n r 



(50) 



holds in light-cone gauge, where we have defined d= \{d — d)- The parton distribution 
depends on x only through the exponential, which allows one to integrate over the x through 
simple partial integrations. The A integration can then be done trivially. So the moments of 
the quark distribution can be expressed in terms of the matrix elements of local operators 
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in the light-cone gauge. However, in this gauge, these gauge-dependent local operators are 
equal to the gauge invariant operators obtained by replacing the partial derivatives in the 
above expression by covariant derivatives, 

g O; in = HO) W {IX1 ■ ■ ■ 7^(0), (51) 

where (• ■ ■) signifies that the indices are symmetrized and the trace has been removed. 
Thus the moments of the quark distribution functions are just matrix elements of the + 
components of the above operators between the initial and final hadron states. We also 
recognize that 

(n • id) q Or + = 2£0+- + • (52) 

This prompts us to define 
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(53) 



After replacing the moments of parton distributions in Eq . (f45"D with matrix elements of 
these operators and interpreting the result as an operator relation, we find the following 
generalized OPE, 

'/^""(IM-f) 

oo n (o n ~ m a ■■■(] \ 

v " / / j / j I I f^2\n J ' ' n—m,m a ■ n,m 

n even=2 m even=0 V \t> / / a=q,g 

oo n / '2 n ~ m q ■ q \ 

+i^^ua/3Q ( ) ^n-m,m a^nfm Mn + ... (54) 

n odd=l m even=0 V \r* ) / a=q,g 

It must be pointed out that the generalized OPE does not have a unique form. One can 
define xb as any dimensionless invariant formed from the external momenta which remains 
finite in the Bjorken limit and expand the amplitude in inverse powers of this variable. This 
will lead to a different set of coefficient functions, but the physical content is the same. The 
choice of which OPE to use is determined by the specifics of the problem at hand. 

Of course, the above expression contains only the contributions to leading order in 1/Q 2 . 
Since the operators O n ^ m are symmetrized and traceless, their rank is n. The mass dimension 
of these operators is just the dimension of the fermion fields plus that of the derivatives, 
or 3 + n — 1. Hence these operators are all twist 2. At the next order in 1/Q 2 one has to 
consider operators of higher twist, which are beyond the scope of this paper. 

VII. SUMMARY AND COMMENTS 

In this paper, we have studied the QCD factorization for deeply virtual Compton scat- 
tering explicitly at one loop and then to all orders in perturbation theory. Our conclusion is 
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that DVCS is factorizable in perturbation theory. This statement has the same level of rigor 
as the ordinary operator production expansion used in deep-inelastic scattering. In fact, 
assuming the generalized OPE with total derivative operators, DVCS can be recovered by 
analytically continue xb variable from xb > 1 region to the point xb = £• The factorization 
theorem guarantees that the Compton amplitudes are finite there, although the one-loop 
calculation indicates that they are not analytic. 

We have also computed the coefficient functions to order a s for the generalized OPE 
including the total derivative operators. For general two photon processes, one has to include 
the longitudinal photon scattering, which has been done in Ref. \T2\, and photon- helictity 
flip amplitude The scale evolution of total derivative operators can best be studied 



using conformally-symmetric operators flT|]. In fact, it has been known for a long time that 
at the leading-log level, the operators of same twist and dimension evolve multiplicatively 
in Gegenbauer polynomial combinations. It is a simple excercise to transform Eq. ([54] ) into 
this basis. 

Note Added: After this work was completed, we learned that the DVCS factorization has 
also been studied by Collins and Freud fl2"4| . Their arguments and conclusions are similar 
to ours. We also learned that some aspects of higher-order corrections to DVCS have been 
considered by Belitsky and Schafer 
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